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Introduction 

T 

X he emergence of imaging genetics to investigate 
the impact of individual genetic variation on brain function 
was presaged by the methodological application of func- 
tional magnetic resonance imaging (fMRI) to schizophre- 
nia research. In 1996, it was first proposed that functional 
magnetic resonance imaging, because of its technical 
advantages over nuclear imaging techniques related to 
enhanced spatial and temporal resolution and noninvasive - 
ness, would enable individual brain phenotype characteri- 
zation for genetic association studies. 1 Further, two seminal 
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reports of gene variation associated with altered brain 
activity served as an initial proof-of-principle and heralded 
the onset of imaging genetics. In 2000, variation in the 
ApoE genotype was reported to be associated with altered 
activity in brain regions affected by Alzheimer's disease 
during a memory task (in hippocampus, parietal, and pre- 
frontal regions) 2 and in 2001, a functional variation in the 
catechol-O-methyltransf erase (COMT) genotype was 
reported to be associated with altered prefrontal activity 
during a working memory task, setting the stage for subse- 
quent investigations of the impact of individual genetic 
variation on brain activity, as detectable by fMRI. 3 
Concurrent with the methodological advances of imaging 
genetics was the conceptual advance of appreciating the 
neuroimaging intermediate phenotype as a manifestation 
of biological risk for a psychiatric syndrome. It was rea- 
soned that susceptibility genes for schizophrenia do not 
directly encode for the clinical syndrome with which they 
are associated, that is, they do not directly encode for psy- 
chopathology, such as first-rank symptoms of hallucina- 
tions or delusions; rather susceptibility genes for schizo- 
phrenia are affecting basic biology — the development of 
brain cells and neural systems that underlie the expression 
of the clinical symptoms of the syndrome. It was posited 
that the effects of susceptibility genes would be more pen- 
etrant, ie, "closer to the gene," at the level of biologically 
based neuroimaging intermediate phenotypes rather than 
at the level of a complex and phenotypically heteroge- 
neous psychiatric syndrome. 4 Neuroimaging intermediate 
phenotypes, akin to cognitive or electrophysiological inter- 
mediate phenotypes, could therefore be used to enhance 
the potential to link genetic variation to a complex psychi- 
atric disorder, such as schizophrenia. 
Hundreds of published articles have ensued, describing 
studies to investigate the association of genetic variation 
with brain activity as pertinent to schizophrenia and 
other CNS disorders. We currently review several critical 
trends in the evolution of the field of imaging genetics 
as applied to schizophrenia research. We then discuss 
where we are poised to go next: innovations in imaging 
analysis and genetics analysis, effective connectivity 
modeling, and polygenic risk models are on the peak of 
the next wave of imaging genetics. 

The neuroimaging intermediate phenotype 

The neuroimaging intermediate phenotype is conceptu- 
ally analogous to an intermediate phenotype for com- 



mon complex medical disorders. It is logical to assume 
that genes would show stronger associations with the 
biological substrates contributing to risk of a disorder, 
with measurable quantitative traits along a pathophysi- 
ologic causal pathway, intermediate to the end complex 
syndrome. Intermediate phenotypes in other realms of 
medicine include lipid level as an intermediate pheno- 
type for heart disease, sodium homeostasis as an inter- 
mediate phenotype for hypertension, and body mass 
index as an intermediate phenotype for diabetes. 5 - 6 We 
favor the term "intermediate phenotype" over the more 
popular term "endophenotype," though the two terms 
are essentially interchangeable. The term "endopheno- 
type" (which was introduced into psychiatric genetics in 
the 1970s 7 ) initially referred to a trait that is "internal" 
that may be discoverable by a "biochemical test or 
microscopic examination," but is not external or overtly 
manifest. Also, the term "endophenotype" does not 
emphasize the concept of intermediacy in pathogenicity. 
Criteria for the establishment of a neuroimaging-based 
intermediate phenotype for schizophrenia, as in other 
fields of medicine, are that the intermediate phenotype: 
(i) is heritable; (ii) is found with increased frequency in 
healthy relatives of ill probands; (hi) exists temporally 
before the onset of the clinical illness in the pathophys- 
iological pathway to the emergence of the clinical syn- 
drome. As expounded in a review by Tan et al, evidence 
for each of these criteria has accumulated for the syn- 
drome of schizophrenia, with cognitive dysfunction often 
integral and assayable at the brain level by task-based 
neuroimaging intermediate phenotypes. 8 Each of these 
criteria are consistent with the assumption that the inter- 
mediate phenotype is genetically and biologically less 
complex than the clinical syndrome, and that genes 
showing association at the syndromal level will show 
greater effect sizes (penetrance) on variation in interme- 
diate phenotypes. 

According to the neuroimaging genetics paradigm, to 
simply demonstrate that a susceptibility gene for schizo- 
phrenia impacts brain function is a necessary but not suf- 
ficient biological proof of a mechanism of susceptibility. 
This is because many, if not most, genes expressed in the 
brain, are apt to have a brain effect of some sort. A sine 
qua non of this proof is to show that the physiological 
intermediate phenotype associated with a susceptibility 
gene for schizophrenia is itself linked to illness risk. To 
make this link, it is necessary to demonstrate that the 
physiological intermediate phenotype is a characteristic 
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of individuals who are at increased genetic risk but do 
not manifest the clinical syndrome. The ideal samples in 
which to demonstrate this are unaffected relatives, eg, co- 
twins, siblings. This has been done for a number of brain- 
associated intermediate phenotypes related to increased 
risk for schizophrenia, including cognitive dysfunctions 
and neuroimaging phenotypes. 914 Thus, the study of 
healthy relatives as a target population is critical for 
establishing the link between genetic association with 
clinical risk, and genetic association with biological risk. 
Having identified a neuroimaging phenotype related to 
increased genetic risk for illness, investigators can ask 
the question of whether genetic variation in a gene of 
interest maps onto the specific phenotype, as an indica- 
tion of its putative neural mechanism of risk. The ques- 
tion arises of which population to choose to conduct this 
test. Neuroimaging studies of only affected subjects is 
confounded by illness-associated epiphenomena that are 
difficult to control, including smoking history, medical 
comorbidities, chronic illness burden, or prolonged neu- 
roleptic exposure. This makes results in patient samples 



difficult to interpret, as the associations may reflect an 
interaction of the gene with any of these epiphenomena. 
Instead, the imaging genetics paradigm to test a specific 
gene-association hypothesis, ie, the association of varia- 
tion in a putative susceptibility gene and brain function 
linked to increased genetic risk, is best performed in 
healthy subjects. Healthy individuals possess common 
at-risk genotypes, but are not themselves symptomatic 
or clinically ill, thereby reducing the effect of confound- 
ing variables. This approach isolates the simple biologic 
effect of the genetic variation on brain function. 

Some neuroimaging caveats 

In contrast to functional neuroimaging, structural imag- 
ing has not yielded robust intermediate phenotypes for 
schizophrenia, although a myriad of structural MRI 
studies have revealed whole brain and region-specific 
volume and other differences between patients with 
schizophrenia and controls. 1517 As discussed elsewhere, 18 
variation in structural MRI measurements may be 



Abnormal behavior reflects abnormal brain function 





Figure 1. The brain imaging intermediate phenotype concept. 
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attributable to artifacts other than the disease itself, 
physiological alterations in brain tissue (eg, tissue per- 
fusion, body fat, or water content), 19 or differences in 
image acquisition and analysis techniques. Alterations 
in body weight, 20 alcohol intake, 21 steroid administration 
and hormonal status, 22 can also change brain volume. 
Of critical importance in evaluating patients with 
chronic schizophrenia, medication can alter brain vol- 
ume, sometimes rapidly; lithium carbonate reportedly 
increased cortical gray matter volume by 3% in 
patients with bipolar affective disorder after 4 weeks of 
treatment, 23 and multiple studies have reported basal 
ganglia volume change with neuroleptics over a dura- 
tion of time intervals, including longitudinal studies of 
first-episode patients. 24 - 25 Furthermore, there has been 
a lack of reliability in reports of the heritability of brain 
structures, with volume reductions in unaffected sib- 
lings inconsistent across brain regions and different 
studies, and not consistently overlapping with genetic 
liability to schizophrenia, decreasing the utility of struc- 
tural volumetric indices as intermediate phenotypes. 26,27 
Increasingly, a growing number of research groups 
have used resting-state fMRI to map brain networks as 
well as diffusion tensor imaging (DTI) to investigate 
white matter abnormalities in patients with schizophre- 
nia compared with controls. These modalities and 
investigations, however, are beyond the current focus 
of imaging genetics in the present review, but have 
been reviewed in depth elsewhere. 28 32 We would 
emphasize two caveats with respect to these 
approaches: (i) the analysis of resting fMRI patterns is 
very sensitive to variation in head motion and to the 
mental state of the subject, making problematic com- 
parisons between ill and well samples 33 ; and (ii) DTI is 
based on highly derived data which are highly suscep- 
tible to many artifacts and biological events unrelated 
to white matter structure, rendering interpretation of 
the results also highly problematic. 34 
In considering activation-based fMRI as applied to 
schizophrenia research there are methodological 
assumptions and limitations to recognize — the same 
limitations of fMRI in general — in its application to 
elucidating other CNS disorders or any cognitive func- 
tion. 35 Firstly, brain mapping is predicated on the 
assumption of a modular organization of the brain, that 
there is a functional segregation with specialized and 
spatially separated modules, with interconnection of 
the entities and functional integration by distributed 



systems. Only according to this assumption can fMRI 
be then employed to reveal hierarchical decomposi- 
tions of brain functional units. Another assumption is 
that a brain structure can be conceptualized as an infor- 
mation processing entity, with an input, a local process- 
ing capacity, and an output. Yet, the traditional cortical 
input-elaboration-output scheme, a correlate of the 
perception-cognition-action tripartite model, may be 
an oversimplification. In contributing to neuronal out- 
put, a change in balance between excitation and inhibi- 
tion may predominate over a more straightforward 
hierarchical connectivity feedforward and feedback 
model. 36 - 37 Further, as for all hemodynamic-based 
modalities, fMRI measures a surrogate signal of brain 
function, but then justifies the assignment of a func- 
tional role to an "active" area, presuming the change in 
signal results from the change of activity of a neuronal 
population. Also, the blood oxygen level-dependent 
(BOLD) signal reflects neuronal mass activity that may 
limit its neurobiological inferences, even while being 
advantageous in some global instances. According to 
one estimate there are about 90 to 
100 000 neurons under 1 mm 2 of cortical surface, 38 or as 
another metric, an unfiltered fMRI voxel contains 5.5 
million neurons, and 2.2-5.5xl0 10 synapses, in contrast 
to traditional microelectrode recordings. Lastly, for 
task-based fMRI neuroimaging studies, using a block 
design, a subtraction method is required to compare a 
task state with an investigator-designed control state, 
requiring a detailed task analysis to determine subtrac- 
tion components and their interactions. Despite the 
above select critiques, fMRI remains the best tool at 
present for gaining insight into brain function, as many 
of these considerations are rendered relatively nonsys- 
tematic because of the principle of every subject serv- 
ing as his/her own control and the comparison of sig- 
nals from one state to another. Thus, fMRI activation 
approaches allow for the ability to test discrete 
hypotheses about definable brain functions, especially 
the impact of genetic variation on these functions. 

Identifying gene associations with 
regions of brain activation 
and neural circuits 

Early studies mapped single genes selected based on a 
measure of candidacy, to an area of brain activation. 
Examples in addition to COMT mentioned above, 
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include the association of a GRM3 variant with ineffi- 
cient prefrontal function during a working memory 
task, 39 or the association of the short (S) variant in a vari- 
able repeat sequence in the promoter of the serotonin 
transporter gene, 5-HTTLPR, with altered amygdala 
activity during an emotionally evocative task. 40 
Illustrative of genetic vulnerability maps of brain func- 
tion based on imaging genetics approaches, Rasetti et al 
reviewed neuroimaging intermediate phenotypes of 
schizophrenia and the gene variants associated with 
them, catalogued by cognitive task. 41 For example, work- 
ing memory tasks have revealed abnormal engagement 
of: dorsolateral prefrontal cortex associated with 
CACNA1C, COMT, GAD1, GRM3, DRD2, KCNH2, 
MTHFR, and RGS4; ventrolateral prefrontal cortex 
associated with COMT, PRODH, and RGS4; thalamus 
and hippocampus associated with NRG1, DAOA/G72, 
and DISCI; and parietal lobe associated with PRODH. 
Cognitive control tasks (designed to challenge executive 
function of goal-directed behavior in the presence of 
conflict) have identified abnormal engagement of the 
anterior cingulate cortex associated with COMT, DRD2, 
and MAO A; of the dorsolateral prefrontal cortex 
(DLPFC) associated with D TNBP1 , DRD2, MAOA, and 
COMT, and of the parietal cortex associated with 
DRD2, and MAOA. 42 Memory encoding tasks recently 
identified abnormal engagement of hippocampus- 
parahippocampus region 43 and association of the hip- 
pocampus with BDNF, COMT, DISCI, GRM3, and 
KCNH2. It is worthwhile to note that most association 
studies of brain function have used single gene variants 
and risk haplotypes emerging from linkage studies and 
more recently genome wide association studies, with dif- 
fering levels of genetic evidence for each candidate gene, 
though there has been no systematic approach to date 
to selecting genes for imaging genetics studies. 
Imaging genetics approaches have progressed to asso- 
ciating gene variants with multiple regions of activa- 
tion, with disease-relevant risk circuits and putative dis- 
tributed functional networks, rather than isolated, 
single regions. After all, brain information processing 
does not occur as discrete activation "blobs," but as 
activity across distributed neural systems and circuits. 
Thus, circuit-based phenotypes would be expected to 
have greater fidelity in showing genetic association at 
the level of brain function, since in principle, the more 
realistic the phenotype, the stronger the genetic associ- 
ation. As schizophrenia is an emergent property of 



neural system function, not isolatable to a singular 
brain region or localized regional defect, but likely 
attributable to network-based neurointegrative deficits, 
neuroimaging and intermediate phenotyping strategies 
have progressed to better understand distributed net- 
works associated with increased genetic risk. To iden- 
tify a functional network or interregional coupling, 
functional connectivity between spatially remote 
regions is inferred based on temporal coherence, by 
identifying regions of coactivation. 44 Statistical analyses 
used for functional connectivity include mapping based 
on seed voxel correlations, principal component analy- 
sis, independent component analysis, and partial least 
squares methods. 

The functional connectivity literature within schizophre- 
nia research has largely focused on PFC connectivity, 
especially the DLPFC and anterior cingulate, and 
DLPFC interaction with the medial temporal lobe, 
specifically the hippocampal formation (FIF), and inter- 
action with the DLPFC-thalamus. 45 For the DLPFC, 
abnormal connectivity has been identified in multiple 
studies in patients with schizophrenia and in high-risk 
subjects 46 " 49 and various genetic associations have been 
established with this putative circuit, during working 
memory tasks. For example during the n-back working 
memory task, COMr-associated activation changes 
mapped onto a prefrontal-parietal-striatal circuit, involv- 
ing the DLPFC, ventrolateral prefrontal cortex 
(VLPFC), the posterior parietal regions, and the stria- 
tum. 50 A PRODH schizophrenia risk haplotype was asso- 
ciated with increased striatal-frontal functional connec- 
tivity, while the protective haplotype was associated with 
decreased striatal-frontal functional connectivity. 51 A 7 
single-nucleotide polymorphism (SNP) haplotype of 
PPP1R1B (encoding DARPP-32) was associated with 
functional coupling and increased activation of the stria- 
tum and prefrontal cortex. 52 An RGS4 variant was found 
to impact frontoparietal and frontotemporal coupling. 53 
A CACNA1C risk SNP was associated with decreased 
prefrontal-hippocampal connectivity. 54 While some func- 
tional connectivity studies have employed a classic inter- 
mediate phenotype strategy, testing for regions of corre- 
lation in affected subjects as well as unaffected relatives, 
many studies appear to query the association of risk sus- 
ceptibility gene variants with correlated regions in 
healthy subjects or affected patients alone, with the puta- 
tive neural circuit then pending validation as an inter- 
mediate phenotype. 
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Psycho-physiological interaction (PPI) analysis is an 
alternative approach to estimating connectivity, and mea- 
sures a regionally specific response in terms of an inter- 
action between a cognitive (or sensorimotor) process and 
activity in another part of the brain. The supposition is 
that the remote region is either the source of afferents 
that confer functional specificity on the target region or 
is activated by efferents that are specifically active during 
the task. PPI, therefore, allows for the exploration of the 
effects of an independent variable (eg, genotype) on task- 
related differences in interregional connectivity. 55 As a 
specific example, combining information about activity 
in the parietal region, mediating attention to a particular 
stimulus, and information about the stimulus, PPI aims to 
identify regions that respond to that stimulus when, and 
only when, activity in the parietal region is high. If such 
an interaction exists, then one might infer that the pari- 
etal area is modulating responses to the stimulus for 
which the area is selective. While this approach offers a 
deeper functional probe of network activity than simple 
time series correlation analyses, it is still based on a cor- 
relation of activity and not on a directional model of 
activity in one region influencing activity in another. 
In one study, during the N back task, DLPFC-HF cou- 
pling was identified in both patients with schizophrenia 
and their unaffected relatives, and associated with a 
ZNF804A risk allele, using seeded connectivity as well 
as PPI approaches. 56 PPI analysis showed a reduction in 
task load-related modulation of coupling between the 
right DLPFC and bilateral HF, in patients and siblings, 
compared with controls. Further, subjects homozygous 
for the risk-associated allele of ZNF804A showed a dis- 
ruption in task-related modulation of right DLPFC-left 
HF coupling in the PPI analysis. The seeded connectivity 
analysis showed similar results to the PPI analysis in 
DLPFC-HF coupling. 

Overall then, functional connectivity analysis offers 
some insight into correlation between different brain 
regions, but is limited in that it does not account for 
directionality, influence, or causality between putatively 
interacting regions; it makes no assumptions about the 
nature of underlying pathways, their structure, nor 
anatomical connectivity. So while correlative methods 
provide a way to characterize neural functional net- 
works by temporal coherence of inter-regional activation 
patterns, it yields neither an understanding of driving 
neural origins nor of the directionality of the observed 
network. 



The next wave of imaging genetics: 
effective connectivity modeling 

In contrast to functional connectivity approaches, effec- 
tive connectivity analyses promise extended insight, refer- 
ring explicitly to the influence that one neuronal system 
exerts over another, and may be used to better explain 
integration within a distributed neural system. Models 
employed in analyzing imaging data to uncover effective 
connectivity are based on regression models, or structural 
equation models, and these models may be linear or non- 
linear. Dynamic causal modeling (DCM) is a type of 
effective connectivity analysis that yields directional, 
pathway information and allows for a quantification of 
the influence of a given neural region over another. 5758 
DCM analysis, introduced in 2003 for fMRI data, is a 
Bayesian framework for inferring hidden neuronal states 
from measurements of brain activity; it is a hypothesis- 
driven approach, requiring an a priori definition of a set 
of interconnected neural areas that mediate a given 
function of interest. 59 DCMs are generative models of 
brain responses, which provide estimates of neurobio- 
logically interpretable quantities including strength of 
synaptic connections among neuronal populations and 
their context-dependent modulation. 60 Causality in 
DCM is based on control theory, ie, causal interactions 
among hidden state variables that are expressed by dif- 
ferential equations that describe how the present state 
of one neuronal population causes dynamics in another 
via synaptic connection, and how these interactions 
change under the influence of external perturbations 
(eg, experimental manipulations) or brain activity. DCM 
tests hypotheses about neuronal mechanisms, allowing 
one to specify a generative model of measured brain 
data, which is a probabilistic mapping from experimen- 
tally controlled manipulations to observed data, via neu- 
ronal dynamics. 

DCM has begun to be applied to imaging genetics. Using 
a DCM approach, distributed circuits that putatively 
underlie working memory — prefrontal-parietal and pre- 
frontal-striatal circuits — were identified in healthy, nor- 
mal subjects, and COMT, DRD2, and AKT1 functional 
variants were associated with the circuits. 61 The goal of 
the study was to engage a hypothesis-driven strategy to 
study component dopamine signaling processes, both D2- 
mediated and non-D2-mediated aspects of dopamine sig- 
naling. The investigators developed a working memory 
task that allowed dissociation of working memory into 
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sub-processes, specifically maintenance of information 
and manipulation of information. In accordance with the 
DCM approach, models of prefrontal-subcortical-parietal 
networks were generated (each model's nodes, connec- 
tions, and inputs were generated) during working mem- 
ory maintenance and manipulation events, and the opti- 
mal model with the highest group Bayes factor was 
determined. The best DCMs for maintenance were pri- 
marily prefrontal-parietal connections, while for manip- 
ulation, the circuit that best fit the data was a prefrontal- 
striatal network. These results fit remarkably well with 
data from nonhuman primates about subprocesses in 
working memory and the principal networks engaged. 
The cortical network engaged during maintenance is pre- 
sumed to be a non-D2 dominated network, and indeed, 
only COMT showed association with activity in this net- 
work. In contrast, the cortical-striatal network is expected 
to be D2-dominated, and all three genes showed effects 
on this network. This study illustrates the greater fidelity 
of genetic association based on more realistic models of 
brain information processing. 

In a study using nonlinear DCM, subjects at high familial 
risk of schizophrenia performed a sentence completion 
task, and the connection strength of the mediodorsal 
(MD) thalamus and inferior frontal gyrus (IFG) was 
investigated, revealing lower connection strength in the 
at-risk subjects. 62 Bayesian Model Selection was used to 
compare the optimal bilinear and nonlinear models, and 
Bayesian Model Averaging was used to assess the con- 
nection strengths with the gating from the MD thalamus 
and the IFG, with nonlinear models providing better 
explanation of the data. In another study, dynamic causal 
models were applied to fMRI data to investigate how 
brain connectivity during an associative emotional learn- 
ing task is affected by different PPP1R1B variants 
(DARPP32-encoding), in healthy subjects. 63 A PPP1R1B 
variant was associated with increased connectivity 
between the inferior frontal gyrus (IFG), amygdala and 
parahippocampal gyrus (PHG), with directionality of the 
connectivity determined to be from the IFG to the PHG. 
In addition to emerging effective connectivity analyses 
by DCM, connectivity is being explored from a more sys- 
tems-level, hierarchical perspective, using graph theory 
metrics to describe the structural and functional compo- 
sition of neural circuits. In graph theory, the correlated 
activity across multiple, distributed preselected brain 
regions can be expressed in terms of a graph, having var- 
ious quantitative parameters, such as nodes, hubs, edges, 



pathway length, and connectivity strength. 64 - 65 The "hubs" 
of these networks correspond to the most highly inter- 
connected neural regions, which typically map to the 
association cortices of the human brain. The strength of 
each node is defined as its average connectivity with all 
other nodes, and the graph's size is defined by the num- 
ber of nodes in the largest connected component; a larger 
graph size indicates fewer disconnected nodes. 6667 
Accumulating evidence suggests that the small-world 
topological properties of brain functional networks are 
altered in patients with schizophrenia. In one study, in 31 
patients with schizophrenia compared with 31 healthy 
controls, functional connectivity between 90 cortical and 
subcortical regions was estimated by partial correlation 
analysis and thresholded to construct a set of unidirected 
graphs. 68 The healthy subjects demonstrated efficient 
small-world properties, whereas topological parameters 
of brain networks — strength and degree of connectiv- 
ity — were decreased in patients with schizophrenia, 
especially in the prefrontal, parietal, and temporal lobes, 
consistent with a hypothesis of dysfunctional integration. 
In another study, in a sample of 203 patients with schiz- 
ophrenia, compared with 259 healthy controls, multi- 
modal network organization was noted to be abnormal, 
as measured by topological and distance metrics of 
anatomical network organization, abstracted from fMRI 
data. 69 Patients with schizophrenia, compared with con- 
trols, demonstrated reduced hierarchy throughout the 
small-world regime, and increased connection distance 
in the multimodal cortical network. The loss of frontal 
hubs and the emergence of nonfrontal hubs was also 
noted, supporting the hypothesis of schizophrenia as a 
dysconnectivity syndrome, impacting the efficiency of a 
frontally dominated hierarchical network of multimodal 
cortical connections. Though the impact of genetic vari- 
ation on network topology based on graph analyses has 
not yet been reported, moderate levels of heritability 
have been found for brain graph topology measured in 
a twin study using EEG, suggesting that genetic variation 
may impact small-world organization and brain graph 
metrics. 70 

The next wave of imaging genetics: 
polygenic risk 

Just as imaging genetics will continue to incorporate 
increasingly sophisticated analytic methodologies, so too 
will imaging genetics evolve to incorporate increasingly 
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sophisticated models of genetic risk, reflective of the 
increasingly apparent polygenic complexity of psychi- 
atric syndromes. Genome-wide association studies 
(GWAS) have indicated a highly significant polygenic 
component of schizophrenia risk, possibly involving up 
to thousands of common alleles of very small effect, at 
the population level. 71 

While early imaging genetics used intermediate pheno- 
types to assess the impact of single gene variants, recent 
studies have increasingly tended towards epistatic mod- 
els of gene interaction. In 2007, Tan et al reported a two- 
way risk variant epistasis: the COMT-Val risk allele, 
associated with reduced prefrontal dopamine, and the 
GRM3 risk allele, related to suboptimal glutamatergic 
function, interacted to give disproportionately inefficient 
DLPFC activation in a working memory task. 72 Further, 
the inefficiency was associated with reduced frontopari- 
etal functional connectivity. Nicodemus et al reported 
the first 3-way interaction using neuroimaging genetics 
to assess the risk susceptibility of the NRG1 molecular 
pathway, finding epistasis between NRG1, and its tyro- 
sine kinase receptor ERBB4, in a 3-way interaction with 
a variant of AKT1P The statistical interaction was bio- 
logically validated by fMRI, in which healthy individuals 
carrying all three at-risk genotypes for NRG1, ERBB4, 
and AKT1 were disproportionately less efficient in 
DLPFC processing than any other combinations of one 
or two at-risk genotypes. Of note, lower-level interac- 
tions were not observed between NRG1, ERBB4, and 
AKT1, suggesting that the interaction, and the NRG1 
pathway, was necessary for the observed fMRI effect of 
inefficiency. Other reports of epistasis in neuroimaging 
genetics include association of variants of with altered 
DLPFC activation, during working memory tasks includ- 
ing DISCl-CIT-NDELl,MTHFR-COMT, and COMT- 
RGS4. 7 " 6 

Imaging genetics is further evolving towards modeling 
increasing genetic complexity, by utilizing a polygenic 
risk score or propensity score of genetic risk for schizo- 
phrenia in fMRI studies. A range of options for con- 
structing a polygenic score may be considered, selection 
of markers according to their P-values in association 
studies, and different methods for weighting markers in 
the score. 77 Only a handful of studies utilizing a poly- 
genic risk score have been reported to date, using both 



functional and structural neuroimaging, and for multiple 
psychiatric syndromes. Walton et al calculated a genetic 
risk score for schizophrenia, the additive effect of 41 
SNPS from 34 putative risk genes, and found a positive 
relationship between the genetic risk score and left 
DLPFC inefficiency during a working memory task. 78 
Holmes et al reported a structural anatomic association 
with polygenic risk for Major Depressive Disorder 
(MDD) in a sample of 1050 healthy young adults with 
no history of psychiatric illness. Using risk scores derived 
from large MDD GWAS analyses, an MDD polygenic 
score was found to be associated with reduced cortical 
thickness in the left medial prefrontal cortex, a structural 
variation that is believed to influence vulnerability to 
MDD. 79 In a third study, increasing polygenic risk allele 
load for bipolar affective disorder (BPAD) was associ- 
ated with increased activation in limbic regions previ- 
ously implicated in BPAD, including the anterior cingu- 
late cortex and amygdala during a verbal fluency task. 80 
So, while a few early imaging genetics studies have 
employed the polygenic risk score, use of the polygenic 
score approach remains to be assessed and validated in 
larger-scale, more robust studies, with an explicit focus 
on schizophrenia and with various models of the risk 
score calculation possible. 

Conclusion 

As highlighted in this review, initial methodological and 
conceptual advances in functional neuroimaging 
launched the epoch of imaging genetics. Since then, 
development of genetic vulnerability maps of the brain, 
identifying neuroimaging intermediate phenotypes of 
schizophrenia and the risk variants associated with 
them, have become a major research industry. While 
imaging genetics to date has led to an increased under- 
standing of schizophrenia pathophysiology and poten- 
tial sites of pharmacologic intervention, a new wave of 
imaging genetics is fueled by even further methodolog- 
ical and conceptual advances. Effective connectivity 
modeling promises to offer causal and directional 
insight into brain networks and circuitry, and polygenic 
risk modeling promises to incorporate genetic models 
reflective of the polygenic complexity of the schizophre- 
nia syndrome. □ 
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Neuroimagenes funcionales y 
esquizofrenia: una mirada hacia el 
modelo de conectividad efectiva y el 
riesgo poligenico 

Se revisan las tendencias mas importantes en image- 
ries y genetica aplicadas a la investigation en esqui- 
zofrenia y se discuten algunas perspectivas a futuro 
en este campo. Gran cantidad de estudios de image- 
nes y genetica han investigado el impacto de la varia- 
tion genetica en la funcion cerebral desde que apa- 
recio el paradigma de un fenotipo intermedio de 
neuroimagenes para la esquizofrenia. Inicialmente 
se postuld que los efectos de los genes susceptibles 
para la esquizofrenia tendn'an una mayor penetra- 
tion a nivel de los fenotipos intermedios de neuroi- 
magenes con base biologica que a nivel de un sin- 
drome psiquiatrico complejo y fenotlpicamente 
heterogeneo. Los resultados de muchos estudios 
apoyan esta hipotesis y la mayon'a de ellos muestra 
variantes geneticas unicas que se asocian con cam- 
bios en la actividad de regiones cerebrates localiza- 
das, como se puede determinar a traves de la selec- 
tion de tareas cognitivas controladas. A partir de 
estos estudios basicos, el analisis de neuroimagenes 
funcionales de los fenotipos intermedios ha progre- 
sado hacia modelos de disfunciones cerebrates mas 
complejos y realistas, incorporando modelos de 
conectividad funcional y efectiva, que incluyen las 
modalidades de interaction psicofisiologica, el 
modelo causal dinamico y las mediciones de la teon'a 
de los grafos. Los enfoques de asociacion genetica 
aplicados a las imagenes y genetica tambien han 
progresado hacia efectos multivariados mas sofisti- 
cados, incluyendo la incorporation de interacciones 
epistaticas de dos o tres vias, y mas recientemente 
modelos de riesgo poligenico. Las imagenes y gene- 
tica constituyen una estrategia unica y poderosa para 
la comprension de los mecanismos neurales de riesgo 
genetico de trastornos complejos del sistema ner- 
vioso central a nivel del cerebro humano. 
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